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Abstract 



Motivated by the observation that there may exist hadronic excitations even in 
the quark-gluon plasma (QGP) phase, we investigate how the properties of quarks, 
especially within the quasi-particle picture, are affected by the coupling with bosonic 
excitations at finite temperature (T), employing Yukawa models with a massive scalar 
(pseudoscalar) and vector (axial- vector) boson of mass m. The quark spectral function 
and the quasi-dispersion relations are calculated at one-loop order. We find that there 
appears a three-peak structure in the quark spectral function with a collective nature 
when T is comparable with m, irrespective of the type of boson considered. Such a 
multi-peak structure was first found in a chiral model yielding scalar composite bosons 
with a decay width. We elucidate the mechanism through which the new quark collec- 
tive excitations are realized in terms of the Landau damping of a quark (an antiquark) 
induced by scattering with the thermally excited boson, which gives rise to mixing and 
hence a level repulsion between a quark (antiquark) and an antiquark-hole (quark-hole) 
in the thermally excited antiquark (quark) distribution. Our results suggest that the 
quarks in the QGP phase can be described within an interesting quasi-particle picture 
with a multi-peak spectral function. Because the models employed here are rather 
generic, our findings may represent a universal phenomenon for fermions coupled to 
a massive bosonic excitation with a vanishing or small width. The relevance of these 
results to other fields of physics, such as neutrino physics, is also briefly discussed. In 
addition, we describe a new aspect of the plasmino excitation obtained in the hard- 
thermal loop approximation. 
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§1. Introduction 



The exploration of phase transitions in quantum chromodynamics (QCD) and the exami- 
nation of the nature of the recently discovered form of QCD matter under extreme conditions, 
such as high temperature (T) and density, are of great interest. One of the central issues in 
this field is revealing the physical content of the so-called quark-gluon plasma (QGP) phase. 
Although it was widely believed for a long time that the QGP phase, even just after the de- 
confinement transition, is composed of weakly interacting quarks and gluons, in accordance 
with the asymptotic freedom of QCD, such a view is now believed to be too naive. In fact, 
it has been being elucidated experimentally and theoretically that the QGP phase near the 
critical point of the chiral and deconfinement phase transitions has a non-trivial structure, 
because the strong coupling nature of QCD plays a significant role. Experiments performed 
at the Relativistic Heavy Ion Collider (RHIC) at BNL are able to probe the QGP phase 
near the critical temperature (Tc) of the deconfinement and the chiral phase transitions. It 
has been argued^^ that the experimental results at RHIC suggest that the created matter 
behaves like a perfect fluid, which imphes that the QGP near Tc is a strongly coupled system, 
and for this reason, it has been called sQGP, where the "s" stands for "strongly coupled" . 
This suggestion seems to be consistent with the results of recent lattice QCD simulations^^ 
showing that the charmonium states survive at higher temperature above Tc of the decon- 
finement transition than originally believed. The existence of hadronic excitations in the 
hght quark sector in the QGP phase has also been suggested, on the basis of the symmetry 
nature of the chiral transition,^)'*) where it was argued and demonstrated using a chiral 
effective model that if the chiral transition is second order or nearly second order, there may 
exist soft modes of the phase transition above Tc that have the same quantum numbers as 
the sigma meson and the pion. 

Quarks and gluons are the basic degrees of freedom in the QGP phase, and therefore 
it is fundamentally important to clarify their properties even apart from possible hadronic 
excitations. In fact, the success of the recombination (or coalescence) model^^ in accounting 
for the quark-number scaling of the so-called V2 parameter of the collective flow seen in 
RHIC experiments strongly suggests that there exist quasi-particles with quark quantum 
numbers in the matter created at RHIC. Therefore it is imperative^) to clarify how the 
strong coupling nature of the sQGP and the quark quasi-particle picture can be compatible. 
The present work is concerned with the properties of quarks in a strong coupling system at 
finite T. 

*^ See also Ref. 5), in which quite different reasoning is given for the possibihty of the color-singlet 
excitations in the QGP phase. 
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It is noteworthy that even in the high-T hmit, where the hard thermal loop (HTL) ap- 
proximation^^ has been established, quarks have some collective nature/°^ as gluons become 
plasmons:^^-*'^^-* The quarks have two branches of spectra, i.e., those of the normal quasi- 
quark and the plasmino,^^) the latter of which does not exist at zero T and density. The 
origin of the plasmino excitation at high T is discussed in Ref. 14),*^ where it is shown that 
the appearance of the plasmino at high T is due to the fact that hole states of thermally 
excited particles and anti-particles can be created. 

In the vicinity of Tc of the chiral phase transition, non-perturbative effects are important, 
and one may expect, on general grounds, that the quarks will possess novel and complicated 
properties. In a previous work,^) the present authors explored the quark spectrum near 
but above of the chiral transition, focusing on the effects of the chiral soft modes. A 
surprising finding in that study was that the quasi-quarks and quasi-antiquarks acquire a 
novel collective nature owing to the coupling of the quarks with the soft modes, and as a 
result, the quark spectral function comes to have a three-peak structure as T approaches Tc7^ 
What causes the three-peak structure in the fermion spectrum? In reality, the soft modes 
correspond to a pronounced peak in the time-like region with a width in the spectral function 
in a bosonic channel: The peak position at the momentum p can be expressed approximately 
as a; ~ ±A/m* (T)^ + with a T-dependent 'mass' ■m'^(T), and as T approaches Tc, m*^{T) 
becomes smaller, together with the width of the peaks. This implies that whenever the multi- 
peak structure of the quark spectral function appears, the soft modes have the character of 
a well-defined elementary boson with a mass, which is found to be comparable with T.^^'**^ 
Thus it is seen that a smaller width of the bosonic excitation is favorable for the multi-peak 
spectral function, and it is feasible that a system composed of a (massless) quark and an 
elementary massive boson, as described by a Yukawa model, may exhibit the multi-peak 
structure in the quark spectral function. 

In this paper, we quantitatively examine the quark properties at finite T using Yukawa 
models with a massive boson coupled to a massless quark. We consider two types of bosons, 
i.e., massive scalar (pseudoscalar) and vector (axial- vector) bosons with a mass m; note that 
because the quark is massless, the scalar (vector) and pseudoscalar (axial-vector) give the 

The origin of a similar spectrum at zero T but high density is discussed in Ref. 15). 
**^ The bosonic excitation corresponding to a peak of the spectral function in the time-like region with a 
small width is a propagating mode with a small damping effect. Such a mode is quite different from those in 
the (color-) super-conducting transition, which are almost diffusive, with their strength concentrated around 
the Fermi surface. Through a coupling with the diffusive soft mode, the quark spectrum can form a 
two-peak structure, -^^^ which causes the pseudogap in the density of states of quarks. -^'^^ Thus it is seen 
that the difference between the numbers of peaks in the quark spectral function comes from the difference 
between the natures of the respective soft modes. 
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same results in the perturbation theory. We calculate the spectral function and the quasi- 
dispcrsion relation of the quark at finite T at one-loop order. We show that a three-peak 
structure in the quark spectral function is formed when T ~ m, irrespective of the boson 
type. We elucidate the mechanism through which the new quark collective excitations are 
realized in terms of the Landau damping of a quark (an antiquark) induced by scattering 
with thermally excited boson, which gives rise to mixing, and hence level repulsion, be- 
tween the quark (antiquark) and the antiquark-hole (quark-hole) in the thermally excited 
quark (antiquark) distribution. Such a mechanism of particle mixing is called "resonant 
scattering". '^'^'-^^^'^'^^ We shall also show that the high-T, weak-coupling limit of the fermion 
spectrum approaches that in the HTL approximation, and thereby show that the plasmino 
excitation obtained in the HTL approximation can be understood as originating from a level 
repulsion between collective quark and anti-quark excitations. 

Our results suggest that the quarks in the QGP phase can be described within an in- 
teresting quasi-particle picture with a multi-peak spectral function, since there may exist 
bosonic excitations in the QGP phase. Furthermore, noting that the models employed here 
are rather generic, it is natural to conjecture that the appearance of the novel three-peak 
structure in the fermion spectrum is a universal phenomenon for fermions coupled to a 
bosonic excitation with a vanishing or small width, irrespective of the type of the bosonic 
excitation. 

This paper is organized as follows. In §2, the spectral function and the (quasi)-dispersion 
relation of the fermion, which are frequently used in the following sections, are summarized. 
We emphasize that the quasi-dispersion relations may not describe the physical excitations 
when the imaginary part of the quark self-energy is large. In §3, we investigate the fermion 
spectrum in a Yukawa model with a massive scalar boson. We show that a three-peak 
structure in the fermion spectral function is formed when T ~ m. We clarify the mechanism 
through which the multi-peak structure in the spectral function is realized in terms of the 
Landau damping and the resonant scattering. The high-T, weak-coupling limit of the fermion 
spectrum is also discussed. In §4, we examine the fermion spectrum in a Yukawa model with 
a massive vector boson. We find that the three-peak structure in the fermion spectral 
function appears for any type of boson. The relevance of the result to QGP physics is briefiy 
discussed. The last section is devoted to a brief summary and concluding remarks. In the 
appendices, some technical details of the calculations given in the text are presented. 



§2. Generalities concerning the queirk spectral function 
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and dispersion relations at finite temperature 



In this section, we summarize general features of the spectral function and the quasi- 
dispersion relations of the Dirac fermion, which we call a 'quark', at finite T. Although the 
contents of this section may be familiar to the reader, this section serves to introduce the 
notation used in the subsequent sections. 

The quark spectral function A{p, to) is expressed as-'^'''^ 

.4(p..) = -lln.G>..).-l g"'^--'-f"'"'--)^° . (2.1) 

where G^{p,uj) is the retarded Green function of the quark, 

with the quark mass m/ and the retarded self-energy a;). The spectral function 

A{pi ui) has the same matrix structure as the Green function. Owing to the lack of Lorentz 
invariance at finite T, A{p, a;) generically has the Dirac structure 

A{p,uj) ^ poip,uj)^^ - py{p,uj)p ■ J + ps{p,uj), (2-3) 

where p = p/\p\ and wc have assumed rotational and parity invariances. The temporal part, 
po{p,uj), represents the quark number excitations of the system. 

When rrif = and TrZ"^ = 0, with Tr denoting the trace over the Dirac index, the system 
possesses chiral symmetry and ps{p,uj) vanishes. In this case, the self-energy is written as 
= i7o7° — i7v7 ■ p, and the Green function appearing in Eq. (2-2) can be decomposed 
using the projection operators A±{p) = (1 ± j^p ■ 7)/2 as 

= G+(p,a;)vl+(p)7° + G4p,u;)A4ph', (2-4) 

with 

G±(p,a;) = iTr[GV^±(p)] = ^. ^| , ^ , (2-5) 

2 u; + irj^\p\- E±{p,u;) 

and E±{p,uj) = {l/2)Tr[E^A±{p)^^] = I^o T The spectral function A{p,uj) can also be 
written as 

= p+(p,a;)vl+(p)7° + p-(p,a;)vl_(p)7°, (2-6) 

where 

P±{P, = --ImGi = po{p, u;) T Pv(p, ^) (2-7) 

TT 
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represents the spectrum of the quark and the anti-quark excitations. Note that E±, G± 
and p± are scalar functions, while E^, and p are matrices in the Dirac space. In a 
non-interacting system, the spectral function of the massless quark is proportional to the 
delta function: p|'^''(p, cj) = 6{lj^ \p\). At zero density, p+(p, a;) and p_{p,u!) have a simple 
symmetric property owing to charge conjugation invariance, 

= -^)- (2-8) 

The poles of G±{p,u!), z = z±{p), represent collective excitations of the quark and the 
antiquark sectors, respectively. They are obtained by solving the equations 

z±{p)T\p\-S±{p,z±{p)) = (2-9) 

in the complex energy plane. The form of the retarded self-energy guarantees that the 
solutions z± lie in the lower-half plane, owing to causality. Equation (2-9) may have several 
solutions for fixed momentum and p±(p, a;) may have a multi-peak structure associated 
with these poles. 

For convenience, we also define the quasi-dispersion relations cu — u!±{p) as the zero of 
the real part of the inverse Green functions: 

^±(P) T IpI = Rer±(p,^±(p)). (2-10) 

Solving Eq. (2-10) is much easier than solving Eq. (2-9) in the complex plane. The functions 
oj±{p) are real numbers, and Eq. (2-10) may have several solutions for fixed p. When the 
condition lml]±{p, u±{p)) <ti is satisfied, and hence the quasi-particle picture in the 

Landau sense is valid, the difference between Rez± and u!± is small. Then u!±{p) can be 
regarded as an approximation of the excitation energy of the quasi-particles and can be used 
to estimate positions of the peaks in p±.*^ We often use the quasi-dispersion relations in this 
paper. 

§3. QuEirk spectrum in a Yukawa model with a massive scaleir (pseudoscaleir) 

boson 

In this section, we investigate the spectrum of a massless quark ip in a Yukawa model 
with a massive scalar (pseudoscalar) boson at finite T. A detailed analysis of a massive 
quark coupled with a massless boson at vanishing momentum is given in Ref. 21). As we 

In this sense, when this condition is satisfied, the peaks of the spectral function and the quasi- 
dispersion relation, which are gauge-dependent in the gauge theories, accurately reflect the position of the 
gauge-independent poles of the Green function. 
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Fig. 1. The quark self-energy at one-loop order. 



find below, the finite boson mass gives rise to unexpectedly interesting complications in the 
quark spectrum, which were not seen in Ref. 21). 

We start from the following Lagrangian composed of a quark and a scalar boson: 

^ = - g4>)^ + ^ - mV') • (3-1) 

Here, g is the coupling constant and m is the boson mass. We note that if the boson is a 
pseudoscalar, g should be replaced with i^^g. As we shall see, however, this replacement 
does not lead to any difference in the results for the quark spectrum. 

3.1. Quark self-energy 
3.1.1. Formulation 

The quark self-energy in the imaginary time formalism at one-loop order, shown in Fig. 1, 
is given by 

/(Pk 
-—^QQ{k,iixin)V{p- k.iujrn- i^n), (3-2) 

where ^o(fc, i^n) = [^<^n7° — ^ • 7]~^ and V{k, ivn) = [{ii^nY' — k^ — m^]~^ are the Matsubara 
Green functions for the free quark and scalar boson, respectively, and a;„ = {2n + l)7rT and 
Vn — 2mrT are the Matsubara frequencies for the fermion and boson, respectively. 

For the coupling with a pseudoscalar boson, a factor of ^75 appears on both sides of Qo 
in Eq. (3-2). However, these factors cancel out because they anti-commute with Qq, and 
thus the self-energy takes the same form as Eq. (3-2). Therefore, the following results and 
arguments in this section hold also for coupling with a pseudoscalar boson, as mentioned 
above. 

In accordance with the standard procedure in the imaginary time formalism, we carry 
out the summation over the Matsubara modes n and the analytic continuation icum ou + ir) 
(see Appendix A for details) . The result is given by 

+ A,ik) -(^^) + /(^/) +A4k)'^^i^^^-J^^^n,^, (3-3) 
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where Ef = \k\, E^ = a/(p — fc)^ + m^, and n{E) = [cxp(£'/T) — l]^"*^ and f{E) = 
[exp{E/T) + 1]^^ are the Bose-Einstein and Fermi-Dirac distribution functions, respectively. 

The self-energy given in Eq. (3-3) has an ultraviolet divergence, which comes from the 
T-independent part, E^{p,uj)t=o = lim^^o ^^(Pi To eliminate this divergence, we 
renormalize this part by imposing the on-shell renormalization condition. (See Appendix C 
for details of the renormalization procedure.) The renormalized self-energy is given by 

with = (a;,p) and — Pfj.p'^. Note that there is no infrared divergence in Eq. (3-3) as 
long as the scalar boson is massive. Since the T-dependent (T 7^ 0) part, a;)T^o = 

IJ^{p,u) — I]^{p,u)t=o, has no divergences, we can evaluate it without regularization (see 
Appendix A). For the numerical calculations of the T ^ part, it is convenient to first derive 
the imaginary part, lmU^{p,u!)T^o, and then determine the real part using the dispersion 
(Kramers-Kronig) relation 

Rer''(p,.),,„ = -ip r ,,-!f]^£!(P4)l££. (3.5) 

where the symbol P denotes the principal value. 

The imaginary part of the self-energy, Imi7^(p, a;), is given by 

d^k 1 



lmU^{p,uj) = -7^9^ J 



(27r)^2E, 

x{A+{k){l + n-f)5iuj-Ef-E,) (I) 



(IV), (3-6) 



+ A4k){n + f)5{u; + Ef-Eb) (II) 
+ A+{k){n + f)6{u;-Ef + Eb) (III) 

+ A_{k){l + n- f)6{u + Ef + E,)Y^ 

with / = f{Ef) and n — n{Eb). Each term in Eq. (3-6), denoted by (I)-(IV), describes a 
definite decay process and has support in the region shown in Fig. 2. The details of the 
physical meaning of each term are discussed in the following. 

Let us denote the external quark as Q and the internal quark as q, as shown in Fig. 1. 
The former is the quasi-quark interacting with scalar bosons, while the latter is a free quark 
in the approximation used here. The physical meaning of each term in Eq. (3-6) is more 
transparent when Eq. (3-6) is rewritten in terms of the transition probabilities. We can 
express defined in §2 as 

U+{p,u;) = u{p)E''{p,tu)u{p), (3-7) 
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Fig. 2. The supports of the terms (I)-(IV) in Eq. (3-6). 



where u{p) is a free quark spinor with the effective mass parameter vP^, which satisfies 
— \/P^)u{p) = with the normahzation u{p)u{p) = 2y/P^. Then, we have 

lmS+(p,Lj) = — / - — ^- — ^ > 27r^ 



X ^[6\p -kf- k,)\M{Q q, b)\'[{l - /)(1 - n) + fn] (I) 

+ [6\p + kf- h)\MiQ, q 6)n/(l + n)+ nil - /)] (II) 
+ [6\p -kj + k,)\MiQ, b ^ q)\'[nil - /) + /(I + n)] (III) 

+ [6\p + kf + h)\M{Q, q, b -> 0)n/n +(!-/)(!- n)]| (IV), (3-8) 

where s is the spin of the quark q and kf = [Ef, k),kb = [Eh, p — k). The amphtudes M are 
given by 

M{Q ^q,b) = M{Q, b^q)= gu{p)u{p), 

M{Q, q^b) = M{Q, q,b^O)= gv{p)u{p), (3-9) 



where v{p) is defined through (|^+ vP^)f (p) = 0, with the normalization v{p)v{p) = — 2vP^. 
The first term in Eq. (3-8) with the statistical factor (1 — /) (1 + ra) describes a decay process 
Q q + b, and that with the factor fn describes the inverse one, q + b —>■ Q. The second 
term with the factor /(l + n) describes Q + q ^ b, and that with the factor n{l — f) describes 
b Q + q. The third and the fourth terms are similar. These processes are schematically 
depicted in Fig. 3. The terms (II) and (III), which involve a thermally excited particle as 
the incident particle, vanish as T — and are known as the Landau damping. They play an 
important role in the structure of the quark spectral function at finite T, as we see below. 
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Fig. 3. The kinetic processes contained in ImI7+(p, oj). (The corresponding inverse ones are not 
shown.) The thick sohd Hne labeled Q represents the quasi-quark, the thin solid lines labeled 
g on-shell free quarks, and the wavy lines labeled h the on-shell scalar boson. The incident 
on-shell particles are thermally excited particles. 



Similarly, it can be shown that the quantity a;) = v{p)E^{p,uj)v{p) contains am- 

plitudes for the processes M{Q q,b), M{Q, q^b), M{Q, q), and M{Q, q,b^O). 
After some manipulations as given in Appendix A, we arrive at the following result 

Imi:±(p,a;) = --^ / dEf[2\p\Ef :f {p' + m'' - u'' + 2uEf)] 



with Ef = (cu^ — — m^)/'2.{uj ± For vanishing momentum, we have ImZ'+ = Imi7_, 



327rp2 

2 

and ImZ"^ is simply given by 

327r ou^ 



x[l + n{u;-Ef)-f{Ef)] 

^ [i\p\Tuj)i7r''T^±uj'')±ujm'']0{p'' -cu''), (3-10) 



Imr«(p = 0,co) = -7°^ ^"^' f^' [n [co,] + / K)] , (3-11) 



with 



= = 2^- ^'-''^ 



3.1.2. Numerical results 

In the upper panel of Fig. 4, we plot lmE-^.{p,uJ)/g^, which is independent of g, for 
vanishing momentum and several temperatures. We see that lmU^{0,uj)/ g'^ vanishes at 
a; = 0, ±m for all T. These values of uj are located in the boundaries of the regions (l)-(IV) 
shown in Fig. 2. The vanishing decay rate at = ±m is due to the suppression of the 
phase space; energy-momentum conservation requires vanishing momentum of the on-shell 
(anti)quark for each process at a; = ±m. Near a; = 0, by contrast, the distribution functions 
in Eq. (3T1) suppress the decay rate, because the on-shell energies of the boson and the 
quark, \u!b\ and la;^!, go to infinity as a; — >• 0. 

From Fig. 4, we see that ImZ'_,_(p, cu) takes finite values only for \uj\/m > 1 at T = 0; 
the corresponding decay processes are represented by (I) and (IV) in Fig. 3. At finite T, 



11 




the Landau damping has a significant effect, and ImZ'+(p, cj) has support in the regions 
(II) and (III) in Fig. 2. The height of the two peaks in these regions increase rapidly as T 
increases, and for T ~ m, the Landau damping comes to dominate the decay rates, owing 
to the processes (I) and (IV) in the energy range shown in Fig. 4. 

Using the dispersion relation given in Eq. (3-5), we can understand the qualitative behav- 
ior of ReZ'_|_(p, uj) from that of the imaginary part: If there is a sharp peak in the imaginary 
part, the real part has a steep increase at the same energy. In order to see this relation, 
we plot Rei7+(p, ct;)/(7^ for p = in the lower panel of Fig. 4. It is seen that there appears 
oscillating behavior around uj = 0, and the amplitude of this oscillation grows rapidly as T 
increases, along with the growth of the peaks in ImZ'+(p, cu). 

The other important property of IJ^{p,uj) shown in Fig. 4 is that both ReZ'+(O,0) and 
ImZ'+(O,0) vanish for all T. The fact that Z''^(0,0) vanishes implies the existence of a 
pole of the quark propagator at the origin (see Eq. (2-9)). (The reason that ImZ'+(O,0) 
vanishes is explained above.) The vanishing of the real part can be easily understood from 
the dispersion relation and the fact that ImZ'+(0,ti;) is an even function of a;. Thus we find 
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Im (p,a))/g 
-0.6 



p/m 




-1.5 -1 



-0.5 0.5 1 1-5 
co/m 



Fig. 5. The imaginary part of the self-energy |Imi7_|_(j>, u;)| for T/m = 1.5. The dashed curves in 
the u-p plane denote u = ±-y/[pp~+n? and u = ±p. 



that ReZ'+(0, u) is an odd function of uj. 

In Fig. 5, we display \lmE^{p,uj)\/g'^ at finite momenta and energies for T/m — 1.5. 
The two separated peaks at p = 0, coming from the Landau damping, tend to overlap, and 
are smeared as p becomes large. 

3.1.3. Asymptotic behavior in the high temperature, weak coupling limit 

Finally, let us consider the asymptotic behavior of S^[p, uj) in the high T, weak coupling 

limit, i.e. gT — * oo and T/m ^ oo, with g'^T 0. 

In this limit, the leading contribution to the imaginary part, lmE^{p,uj), comes from 

the terms which explicitly include the factor in Eq. (3-10). Thus, the asymptotic form of 

ImZ'^(p, a;) is given by 



Imr^(p,a;)|T- 



32 



ujp ■ 7 

p2 



(3-13) 



The real part in this limit is easily calculated using Eq. (3-13) and the dispersion relation 
given in Eq. (3-5). We then obtain the following asymptotic form of E^{p,uj): 



mi^Qlp, Lu) + mi'^T-^ ( 1 - r^Qip, uj] 
\P\ \P\ \ \P\ 



with the thermal mass my = gT/A and 



\P\ 



UJ — \p\ 



i-eip'-u^ 



(3-14) 



(3-15) 
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Equation (3-14) is the same as the fermion self-energy of QED and QCD in the HTL approx- 
imation/°^ apart from the definition of the thermal mass. Therefore, the quark spectrum 
in the Yukawa model in the limit gT oo [T/m — >■ oo) with — >■ is equivalent to the 
well-known spectrum in the gauge theories. ^^-^ 

We note that the analytic structure of the self-energy is qualitatively different from that 
at intermediate T, i.e. for T ~ m; there is no pole of the propagator at the origin in 
this limit, because ReZ'+(0,c<j — * 0) behaves as ~ 1/uj and hence never vanishes, while 
ImZ'+(0,u; — > 0) vanishes. This change in the self-energy is reffected in that of the peak 
structure of the spectral function, as discussed in the next subsection. 

3.2. Quark spectral function at various temperatures 

In this subsection, we examine how the quasi-particle picture of the quark changes at 
finite T by studying the quark spectral function p±{p,uj) and the quasi-dispersion relation 
oj±{p). We fix the coupfing constant as 51 = 1 throughout this subsection to elucidate the T 
dependence of the spectrum. The g dependence will be discussed in a subsequent subsection. 
Because the only dimensional parameter in Eq. (3-1) is the boson mass m, we scale all the 
dimensional parameters T, uj and p by m. 

3.2.1. Quark spectral function at zero temperature and the hard thermal loop approximation 

Before presenting the numerical results, we briefiy review the quark spectrum at T = 
and high T fimit in our model. The analytic form of the quark self-energy at T = is given 
in Eq. (3-4). In this case, the poles of the quark propagator are on the light cone, a; = ±|p|, 
because we have employed the on-shell renormalization condition to derive the self-energy 
appearing in Eq. (3-4). The quark spectral functions p±(p, tf) are then given by 

= Z(0) • 5{u T \P\) + pf'''\pM. (3-16) 

where the T-dependent residue Z{T) is exactly unity at T = under the on-shell renor- 
malization condition, and the continuum part takes finite values for |a;| > ^Jp^ + rri^. 

In the high T, weak coupling limit, the quark spectrum approaches that calculated in the 
HTL approximation in the gauge theories with the thermal mass tut = gT/A, as does the 
quark self-energy (see the previous subsection). In this limit, both a;) and p_ (p, u) have 
two delta-functions corresponding to the normal quasi-particle and plasmino excitations, in 
addition to the continuum part in the space-like region^^^'^^^ (see Eqs. (3-14) and (3-15)). In 

*^ Some readers may be uneasy about the violation of the sum rule for the spectral function, 
dLup±{p, w) = 1. It should be noted here that in relativistic field theory, this sum rule does not necessarily 
hold after the renormalization. 
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T/m = 0.4 



T/m = 0.8 
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0.4 -0.8 
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Fig. 6. The quark spectral function p^{p,uj) for T/m = 0.4 and T/m = 0.8. 



the high T hmit with fixed g, these two peaks have widths of the order of grriT^^^ 




T/m=0 
T/m=0.4 
T/m=0.6 
T/m=0.8 



0.1 0.2 0.3 0.4 0.5 

p/m 

Fig. 7. The quasi-dispersion relation ijo+{p) for several values of T. 



3.2.2. Quark spectral function at low temperature 

In Figs. 6 and 7, we plot the spectral function p_|_(p, u;) and the quasi-dispersion relation 
oj+{p) of the quark part for relatively low temperatures, i.e., T < m. The anti-quark parts 
p_ and uj^{p) can be determined from the symmetric properties given in Eq. (2-8) and 
the relation uj^{p) = —uj+{p). At T/m = 0.4, the spectral function p_|_(p, cu) and uj^{p) 
are almost the same as those at T = 0; p+{p,u!) has sharp quasi-particle peaks around 
uj = \p\, and uj+{p) ^ At T/m = 0.8, the quasi-particle peaks are clearly deformed 
near \p\/m = 0.15. The quasi-dispersion relation also deviates from the free one near 
\p\/m = 0.15. We note that this relation exhibits an unphysical acausal behavior; the 



15 



T/m = 0.8 



T/m = 1.4 




imaginary part of the self-energy is large in this momentum region, as shown in Fig. 8, and 
thus the quasi- dispersion relation becomes unphysical in this region. 

To understand the behavior of p+(p, to) and uj+{p) depicted in Figs. 6 and 7, we show the 
real and imaginary parts of E^{p,uj) at T/m = 0.8 and several momenta in the left panel 
of Fig. 8. We see that there are two peaks in ImZ'+(0,ti;) for positive and negative energies, 
corresponding to the terms (II) and (III) in Eq. (3-6), respectively. It is seen that ReZ'+(p, u) 
exhibits a steep rise in the two regions corresponding to these peaks. We also plot the lines 
u — \p\, i.e. the l.h.s. of Eq. (2-10), by the dotted lines in Fig. 8. Since the quasi- dispersion 
uj+{p) is given by Eq. (2-10), uj^{p) corresponds to the points of intersection of the dotted 
lines and ReZ'+(p, tu). We see that the crossing points depend strongly on the shape of 
ReZ'+(p, cj). Near u/m = 0.1, the rapid change of ReZ'+(p, cj) causes the quasi-dispersion 
relation to deviate from the free quark dispersion, as shown in Fig. 7. 

3.2.3. Quark spectral function at intermediate temperature 

As T is raised further so that T becomes comparable with m, p±{p, uj) and uj±{p) change 
drastically; note that the bosons become thermally excited with a considerable probability at 
such a temperature. In the left panels of Fig. 9, we plot p+(p, uj) for T/m = 1.0, 1.4 and 1.8. 
We see that for T/m = 1.0, the quasi-particle peak starts to split near \p\/m = 0.15. Also 
in this vicinity, there appear broad peaks both in the positive and negative energy regions 
for lower momenta. Thus the spectral function in the low-momentum region has three peaks. 

A remark is in order here. Among the three peaks, that near the origin has no width, 
and thus it is difficult to read off the strength from Fig. 9. Here we note that the strength at 
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T/m = 1 




T/m = 1.4 




T/m = 1.8 




Fig. 9. The quark spectral function p+{p, co) (left) and the quasi-dispersion relations uj±{p) (right) 
for T/m = 1.0, 1.4, and 1.8 from the top. In the right panels, the (u;,p) dependence of p±{p,uj) 
is represented by the color density. 0;+ (aj_) and p+ (p-) are shown in the right (left) halves 
of these panels. Note the direction of the momentum scale in the left half-plane is opposite to 
that of the right half-plane. 
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Fig. 10. The temperature dependence of the residue of the pole at the origin {uj = \p\ = 0). 

the origin is actually given by the residue of the Green functions, G±: The residue is simply 
given by 

Z{T) = [1 - dReE±/du;]-\ (3-17) 

since the imaginary part, ImZ'-|-(p, u;), vanish at the origin. The T dependence of Z{T) is 
shown in Fig. 10. It is seen that Z{0) = 1, as mentioned previously, and Z(T) decreases 
gradually as T increases but is still considerably large for T/m ~ 1, and hence the three-peak 
structure of the spectral function is realized in this temperature region, as mentioned above. 

In the high T limit, Z(T) vanishes, and only two peaks remain in p±{p,uj). As we 
see below, these two peaks correspond to the normal quark quasi-particle and plasmino 
excitations obtained in the HTL approximation. 

In the right panels of Fig. 9, we plot the quasi-dispersion relations and the {uj,p) de- 
pendence of the spectral function; the latter is represented by the (color) density. In these 
panels, cj_|_ and p+ are shown in the right (left) halves of the figures. Note that 

the direction of the momentum scale in the left half plane is opposite to that in the right 
half plane. At T/m = 1. there is one quasi-dispersion curve for the entire momentum region. 
As T increases, the quasi-dispersion curve bends greatly forming 'back-bending' parts, and 
there eventually appear multi- valued quasi-dispersion relations for some momentum regions. 
At T/m = 1.4, the quasi-dispersion relation is five- valued at low momenta; note that this 
back-bending quasi-dispersion relation is acausal and unphysical. Only three of these values, 
however, are accompanied by peaks of the spectral function, as the quasi-dispersion relation 
a;+(p) near the back-bending region does not form a peak in p+(p, a;) and is unphysical. The 
back-bending feature of the quasi-dispersion relation becomes more prominent as T increases 
further, as shown in the bottom-right panel in Fig. 9. 

To understand the peculiar behavior of p+(p, a;) and u!+{p) displayed in Fig. 9, we show 
Z'_|_(p, cu) at the intermediate temperature T/m — 1.4 in the right panel of Fig. 8. We see 
that the two peaks in lmE^{p,u) become sharper than those for T/m — 0.8, shown in 
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the left panel, and the oscillatory behavior of the real part for ~ also becomes more 
prominent. In order to determine the quasi-dispersion relation uj_^.{p), i.e., the solutions of 
Eq. (2-10), from the figure, we plot the line u — \p\, i.e. the l.h.s. of Eq. (2-10). In the top 
panel, we see that there appear five crossing points. The crossing points with the second 
and fourth largest u, however, are located at the energies of the peak of |ImZ'+(p, cj)!, and 
hence the spectral function does not form a peak there. For large momenta, the number of 
crossing points decreases, and eventually only one crossing point remains. 

3.2.4. Quark spectral function at high temperature 

T/m = 2 T/m = 3 T/m = 5 




1 0.5 0.5 1 1 0.5 0.5 1 1 0.5 0.5 1 

p/mj p/mj p/rrij 



Fig. 12. The quasi-dispersion relations uj±{p) for T/m = 2,3 and 5. All the values are in units of 
the thermal mass, mx = gT/A. The dashed line represents the dispersion relation in the HTL 
approximation . 

Here we examine the quark spectrum at higher temperatures. We plot p_|_(p, cj) and 
'^±{p) foi' T/m = 2.5 in Fig. 11 as an example of the high-temperature case, i.e. T ^ m. 
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We see that two well-formed peaks appear in the positive and negative energy regions, while 
the peak near the origin immediately disappears as the momentum becomes finite. The 
curve representing the quasi-dispersion relation in the right panel is quite complicated. It 
should be noted that, however, the back-bending behavior of the quasi-dispersion relation 
is unphysical, and does not correspond to a peak of the spectral function, as mentioned 
previously. We plot the quasi-dispersion relation at several temperatures, T/m = 2, 3 and 
5, in the high-T region in Fig. 12, together with that in the HTL approximation with the 
thermal mass rriT = gT/A. From the figures, it is seen that our quasi-dispersion relation 
approaches that obtained in the HTL approximation as T increases; the two curves of uj^{k) 
approach the quark quasi-particle and plasmino dispersion relations, respectively, with the 
thermal mass mj- = gT/A.*^ 

Here we note that there appears a continuous spectral bump near c<j = in the space- 
like region in the HTL approximation, which is reminiscent of the spectral peak seen at 
intermediate temperatures in the space-like region near uj = 0. We stress, however, that the 
physical origins of these peaks are completely different. The peak appearing at intermediate 
temperatures is due to the pole of the Green function, as shown in the previous subsection, 
while the bump in the HTL approximation is due to the Landau damping, and the quark 
propagator calculated in this approximation does not have a pole at a; = 0. 

3.3. Discussion of the three-peak structure of the spectral function for T m 

We have seen that the appearance of the three-peak structure for values of T comparable 
with the boson mass m is the most characteristic feature in the quark spectral function 
caused by the coupling with the massive boson. In this subsection, we attempt to elucidate 
the physical origin of the multi-peak structure. This discussion proceeds analogously to 
that^^ given for the origin of the three-peak structure caused by the soft mode of the chiral 
phase transition. 

To understand the physical mechanism responsible for the three-peak structure of a;), 
we first recall that there develop two peaks in ImZ'+ at a positive and negative energy in 
this temperature region; the peak height here increases with T. As discussed in §3.1, these 
two peaks correspond to the decay processes depicted in (II) and (HI) of Fig. 3, which are 
both Landau damping. The process (II) is the annihilation process of the incoming quark Q 
and the thermally excited antiquark into a boson in the thermal bath, Q + q ^ b, and its 
inverse process. Two remarks are in order here. First, the disappearance of an anti-quark 
implies the creation of a 'hole' in the thermally excited anti-quark distribution. Second, 

In fact, the spectrum in Fig. 12 does not exactly coincide with that in the HTL approximation because 
the couphng is not small. However, it is known that their difference is small at one-loop order. 
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the creation of bosons in a thermal bath is enhanced in comparison with the case in vacuum 
by a statistical factor of 1 + n, which becomes large when T is comparable to m. Thus, we 
see that the process (II) causes a virtual mixing between the quark and 'anti-quark hole' 
states through the coupling with the boson in a thermal bath, *^ and as a result, the mixing 
is enhanced when T/m ~ 1. 

The process (III) is another decay process of a quasi-quark state Q, which is now un- 
derstood to be a mixed state of quarks and antiquark-holes, into an on-shell quark via a 
collision with a thermally excited boson: Q + b q and its inverse process, q ^ Q + b. As 
is seen from the left panel of Fig. 2, the energy involved in probable processes occurring at 
small momenta is negative for the process (III). These processes again give rise to a mixing 
of a quasi-quark and an anti-quark hole state. As seen from the Bose-Einstein distribution 
function n, thermally excited bosons are abundant when T approaches m. 

A similar interpretation applies to the anti-quark sector, ImZ'_(p, cj), if the quark and 
anti-quark hole are replaced by an anti-quark and a quark hole, respectively. Thus, the 
process corresponding to (II) in Fig. 3 is Q + q ^ b, where q denotes a thermally excited 
quark, for instance. For this reason, from this point we only consider the quark sector. 

The mechanism for the mixing of the quark and hole state can also be characterized 
as a resonant scattering ^^'^'> which was originally introduced to understand the non-Fermi 
liquid behavior of fermions just above the critical temperature of the superconducting tran- 
sition. ^^^'^'^^ In fact, we have seen that the process (II) in Fig. 3 includes a scattering process 
of the quark by a massive boson, thus creating a hole state in the thermally distributed 
anti-quark states: Q qh + b. Such a process is called resonant scattering. ^'^^ Note that 
the most probable process for finite T involves the lowest energy state of the boson, i.e., a 
rest boson with a energy m. The energy conservation law in the most probable case for the 
above process is ujqip) + u!g{—p) = m, or equivalently, 

ujq{p) = m - ujg{-p). (3-18) 

This equation actually represents the energy-momentum relation for the particles involved 
in the state mixing. Thus we see that the physical energy spectrum is obtained as a result 
of the level repulsion between the energies u!q{p) — \p\ and m — u!q{—p) = m — |p| in the 
perturbative picture. This situation is depicted in the upper-right part of the right panel in 
Fig. 13. 

Similarly, the process (III) in Fig. 3 includes the process g — > Q + b, and the energy- 

We remark that quark number conservation is not violated with this mixing, because an 'anti-quark 
hole' has a positive quark number. 
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Fig. 13. Typical peak position of the spectral functions in the case of level mixing at the origin 
(left) and away from the origin (right). The long dashed curves show that the strength of 
the spectrum is becoming weaker. The solid line represents the free quark and antiquark 
dispersion relations, and the dotted line the antiquark hole and quark hole dispersion relations. 
Level repulsion takes place at the intersection point of these two lines. 

momentum conservation law for this process for the most probable case is 



Thus, the physical energy spectrum is obtained as a result of the level repulsion between the 
energies —\p\ and —m + \p\. This situation is also depicted in the lower-right part of the 
right panel in Fig. 13. 

We thus find that at temperatures satisfying T/m ~ 1, owing to the finite boson mass, 
the level repulsions occur far from the origin, and then the quasi-dispersion relations are bent 
twice, or two gap-like structures in the quark spectrum are formed at positive and negative 
energies, as shown in the right panel of Fig. 13. 

It is interesting to consider the high temperature limit, T 3> m, or m/T ~ 0. In this 
case, the effect of the boson mass can be ignored, and the resonant scattering occurs only 
once at the origin (a; = |p| = 0), since the energy levels which are to be repelled cross only 
there. Then the situation becomes that represented in the left panel of Fig. 13. 

It has been shown''^ that the quark spectral function possesses a three-peak structure 
near Tc of the chiral transition when the chiral soft mode"^-* is incorporated into the quark 
self-energy. As was described in the Introduction, the soft modes behave like a massive 
elementary boson with a mass m*(T) as T approaches i.e. Wgoft ~ ^/p"^ + ^*a{TY, and 
hence the quark spectra studied in Ref. 7) are essentially the same as that treated in this 
section. We also note that as T is lowered toward m*^{T) tends to vanish, and hence the 
ratio T /m%{T) becomes large. Thus the quark spectrum approaches that in the T/m —> oo 
limit of the present work.^^ 




(3-19) 
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9 = 0.3 , T/m = 2.8 



g = 2,T/m=l 




Fig. 14. The quark spectral functions p+(p,ti;) for g = 0.3, T/m = 2.8 and g = 2, T/m = 1. 



3.4. Quark spectral function for various coupling constants 

We have derived the quark spectrum with the fixed couphng g = 1, and we have found 
novel three-peak structure at intermediate temperatures. How do these features change as 
g is varied? In the previous subsections, we pointed out that the Landau damping, which 
induces the two peaks in ImZ'^(p, a;), plays an essential role in creating the three-peak 
structure in the quark spectrum. Because S^{p, u) / g"^ is independent of g at one-loop order, 
we conjecture that the three-peak structure in the quark spectrum would not be affected 
qualitatively even if we varied the coupling g. We confirmed numerically that this is indeed 
the case: There appears similar three-peak structure in the quark spectrum for all values 
of g that we considered, although the temperature at which the clear three-peak structure 
appears depends on g. 

In Fig. 14, we plot the quark spectral function p+(p, a;) with g = 0.3 and 2, for values of 
T where a clear three-peak structure is seen. The spectral function possesses a single-peak 
structure at values of T lower than in the case considered in Fig. 14, and it approaches the 
two-peak structure at higher T. It is notable that for larger g, the three-peak structure 
appears at lower T and the distance between the peaks becomes larger. 

§4. Quark spectrum in a Yukawa model with a massive vector (cLxial-vector) 

boson 

In this section, we investigate the quark spectrum using a Yukawa model with a massive 
vector (axial- vector) boson and show that the three-peak structure of the quark spectral 
function is obtained. 

We start from the following Lagrangian, composed of a massless quark ip and vector 
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boson field V^: 



£ = V5(^^-i^7%)V + >Cv. (4-1) 



We consider the following simplest form for the Lagrangian of the massive vector field, C^: 

£v = -\f^.F^'' + \m^V,V^, (4-2) 

with the field strength F^,^ = — d^V^j,. 

As in the case for the scalar boson, we find that the replacement of the vector boson 
with an axial-vector boson does not lead to any difference in the quark spectrum in the 
perturbation theory in this work. 

4.1. Quark self-energy 

At one-loop order, the quark self-energy in the imaginary time formalism is given by 

E{p,iuj^) = -g'^T22 J j^-j^l^Go{k,iUn)Y'^f^AP- ^^i'^m- ii^n): (4-3) 



with the Matsubara propagator for the massive vector boson, 

PfiP'^ — 

where — {ivn^p)- The propagator for this vector boson, which is obtained from Eq. (4-2), 
is called the Proca propagator. It is known that the massless hmit of Eq. (4-4) cannot be 
taken, because of the lack of gauge invariance of Eq. (4-2). As we see below, this property 
causes a problem in the the high temperature limit, T/m — > oo, because this limit also rep- 
resents the massless limit with fixed T. One should describe a boson field in the Stuckelberg 
formalism,^^) which has a massless limit, if we study the spectrum at high T. In this work, 
because we explore the qualitative effect of the phenomenological massive vector boson on 
the quark spectrum at intermediate T, we do not consider this problem associated with the 
Proca propagator. 

For couphng with an axial- vector boson, the self-energy has the same form as Eq. (4-3), 
because the 75 matrices in the vertices cancel out in the case of a massless quark propagator 
Qq. Therefore, the following results hold also for the coupling with an axial- vector boson, as 
mentioned above. 

Carrying out the summation over the Matsubara modes in Eq. (4-3) and applying the 
analj^ic continuation iuj^ cu -\-ir) (see Appendix B for details), we obtain 

E {p,u;} - g 2^JJ ^ — [^g,. - ^ ^ _ - tE,' 

(4-5) 
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with = {uo — Ef.p — k). In order to eliminate the uhraviolet divergence in Eq. (4-5), 
we employ the same strategy as in the previous section: We first divide Eq. (4-5) into T- 
independent and T-dependent parts as E^(j),uj) = E^{p,uj)t=o + ^^{p,'-^)Ty^Q and calculate 
E^{P,uj)t=o while imposing the renormalization condition. 

In the present case, however, unlike in the previous section, the T-independent part, 
E^{p,uj)t=o, includes a divergence in the term proportional to and the on-shell renor- 
malization condition alone cannot eliminate the divergence (see Appendix C). This is due 
to the bad ultraviolet behavior coming from the Proca propagator Dj^^^. Here we impose the 
renormalization condition d^E^{p)/d;^^\^o = to eliminate this divergence and obtain the 
following renormalized form: 



-—-2 2—^ 
6 P^ / 



,P2 + 2m2 (P2 - ^2^2 



m 



(4-6) 



The T-dependent part E (p,uj)t^o can be calculated using the same procedure as in the 
previous section; i.e., we first calculate the imaginary part ImZ'^(p, a;)r^o and then derive 
the real part using the dispersion relation, Eq. (3-5). 
Prom Eq. (4-5), we obtain 

Imr«(p,w) = -TT^ / — — [g^^ 



x{-(l + n - f)5{u -Ef- Eh) -{n + f)5{u + Ef - E,) 

+(n + f)5{uj -Ef + Eb) + {l + n- f)S{u; + Ef + Ef,)}, (4-7) 

with n = n{Ei,) and / = f{Ef). The four terms in Eq. (4-7) have the same statistical factors 
and delta functions as those in Eq. (3-6), which means that the decay processes of the quasi- 
particles described by these terms can be understood diagrammatically from Pig. 3. The 
region in the energy-momentum plane where each term has a finite value is the same as in 
the case of Pig. 2. 

Por p = 0, we obtain 

Imr±(p = 0,a;) = ^ — coth — — + tanh 



327r m2 u;3 V ^Tuj ATu; 

(4-8) 

Equation (4-8) differs from Eq. (3-11) in the Yukawa model with a scalar boson by a factor 
of (a;^ -|- 2m^)/m^. Por small energy, i.e. a; < m, this term approximately gives an overall 
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factor of 2. Therefore, we expect that the two-peak structure of lmE±{p = 0,u!) caused by 
the Landau damping at < m and, hence, the three-peak structure in the quark spectral 
function is not quahtatively altered even in the present case. For large u, however, the factor 
behaves as tu^/m^, which makes Eq. (4-8) much larger than Eq. (3-11) for \uj\ ^ m. This 
leads to a difference between the scalar and the vector bosons at very high T. 

4.2. Quark spectrum at intermediate temperatures 

In Fig. 15, we plot the imaginary and real parts of uj) for p = at T/m = 0, 0.5, 1 

and 1.5. We fix the coupling constant at g = 1. As can be deduced from the above discussion, 
the qualitative features of both parts for |c(j|/m < 1 are quite similar to those in Fig. 4: There 
are two clear peaks in ImZ'+ for \uj\/m < 1 and oscillating behavior in ReZ'+(p = 0,u!), and 
they grow rapidly as T increases. For \uj\/m > 1, it is seen that |ImZ'_,.(p = O,^^)! grows 
more rapidly than in the Yukawa model with the scalar boson shown in Fig. 4. 

In Fig. 16, we plot |Imi7_|_(p, in the energy-momentum plane for T/m = 1.5. We see 
that the qualitative features of |ImZ'_|_(p, c<j)| are again almost the same as those in Fig. 5 
for \u!\/m < 1. 
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(o/m 

Fig. 16. The imaginary part of the self-energy, \lm.U^{p,Lo)\, for T/m = 1.5. 

In Fig. 17, we plot the quark spectral function for T/m — 0.4,0.8,1.0 and 1.2. It is 
seen that there appears a three-peak structure at intermediate temperatures, as found in the 
previous section. This is quite natural, and it would be expected from the behavior of the 
self-energy. A clear three-peak structure is formed at lower values of T than in the scalar 
boson case. This is because the term (a;^ + 2m^)/m^ in Eq. (4-8), which gives approximately 
a factor of 2 for a; < m, tends to enhance the coupling constant g. As discussed in §3.4, the 
temperature at which the clear three-peak structure appears decreases for larger g. 

4.3. Quark spectrum at high temperatures 

Next, we plot the quark spectral function at higher temperatures, T/m = 5 and 10, 
in Fig. 18. We find that the energy at the quasi-particle peak converges to a value near 
uj/m = 2.5 as T/m increases for p = 0. It is also noteworthy that the two peaks with finite 
thermal mass become obscure as T increases. 

In order to understand this behavior of p+(p, u;) at higher T, in Fig. 19 we plot the self- 
energy for T/m — 0, 5, 10, and 15. In the right figure, we plot the line Rei7_|_ = a; in order 
to see the quasi-dispersion relation. At T = 0, there exist three quasi-dispersion relations at 
a; = and near u = ±10m. The latter two solutions are unphysical.*-* In Fig. 19, it is also 
noteworthy that ReZ'_|_ changes sign near ou — 2.5m for higher T. In fact, we can analytically 
prove that the zero of ReZ'4. approaches uj — \f^m in the high T limit (see Appendix D). 
Due to this property, the solution of the quasi-dispersion relation corresponding to the two 
peaks in Fig. 17 approaches — \f^m at p = in the high T limit. 

*^ These quasi-dispersion relations do not form a peak in the quark spectral function, because jlmZ'^l 
takes a large value there. We have confirmed that they also exist for very large values of a; in the model 
with a scalar boson in §3. 
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Fig. 18. The quark spectral function p+(p, w) for T/m = 5 and 10. 



§5. Summary and concluding remarks 

We have investigated the quasi-particle picture of a fermion, which we call a "quark" in 
this paper, employing Yukawa models composed of a massless quark and a massive boson 
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Fig. 19. The imaginary and real parts of the self-energy, X!^{p,uj), at zero momentum for T/m 
5, 10 and 15. 



at finite temperature (T). We have considered four types of massive bosons, i.e., scalar, 
pseudoscalar, vector and axial-vector bosons, with a mass m, though there is no difference 
in the quark spectra for the scalar (vector) and pseudoscalar (axial-vector) bosons, because 
the quark is massless. The Proca formalism was employed for the vector boson. For both 
models, the quark self-energy was computed up to one-loop order, and the spectral function 
and quasi-dispersion relation of the quark were studied. We have found that the quark 
spectral function at low frequency and low momentum is significantly different from the 
free quark spectrum, and it has a three-peak structure at intermediate values of T, i.e., for 
T/m ~ 1. We have also found that the number of branches of the quasi-dispersion relation 
does not necessarily coincide with that of the peaks in the fermion spectral function when 
the imaginary part of the quark self-energy is large. We have shown that the two of the 
three peaks have finite thermal masses and approach the normal quasi-quark and plasmino 
excitations in the HTL approximation in the high T limit, while the strength of the peak at 
the origin becomes weaker and disappears in this limit. In the Yukawa model with a vector 
boson, we have found that the behavior of the spectral function is the same as that in the 
Yukawa model with a scalar boson up to T m. In the high T limit, the spectrum in the 
Proca formalism does not approach that in the HTL approximation. 

We have discussed the fact that the three-peak structure originates from the Landau 
damping, that is, the scattering processes of a quark and an antiquark hole of a thermally 
excited antiquark, and of an antiquark and a quark hole of a thermally excited quark. 
They cause the formation of energy gaps in the quark spectrum, owing to the level mixing 
between the quark (antiquark) and the hole of the thermally excited antiquarks (quarks). 
This leads to a level repulsion or gap in the resultant physical spectrum. These mixings can 
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be understood in terms of the resonant scattering^^-*'^"-* of the quasi-quarks off the bosons. In 
particular, owing to the mass of the boson, the level repulsion due to the resonant scattering 
occurs twice at different points in the energy-momentum plane, leading to the three-peak 
structure of the spectral function. This contrasts with the case of the quark spectrum coupled 
to a massless boson, such as a gauge boson, in which the resonant scattering occurs only 
at the origin in the energy-momentum plane and then leads to only two peaks, as in the 
spectrum in the HTL approximation. 

Since there may exist bosonic excitations (even apart from gluons) coupled to quarks 
in the QGP phase,^-*'^^-*'^^-* it is plausible that the quark spectrum would have a three- 
peak structure in the QGP phase if the boson has a "mass" of the order of T; if there 
exist several bosonic excitations with "masses" close to T, the quark spectral function may 
become a superposition of three-peak structures with various heights. Although the quark 
spectral function in QCD is a gauge-dependent quantity, the complex pole of the quark 
Green function is gauge invariant. For this reason, wc conjecture that the peak structure of 
the quark spectral function reflects the complex pole of the Green function and will be seen 
in gauge theories, for example, QCD. 

One interesting example of bosonic excitations which may affect the quasi-particle picture 
of the quarks in the QGP phase is the soft modes associated with the cliiral transition.'^' 
The chiral soft mode is unique, because its effective mass m*(T) may become as small as a 
few hundred MeV, and thus it can eventually be of the same order as T near of the chiral 
transition if the phase transition is second order or weakly flrst order. In fact, it has been 
shown^' that the chiral soft modes give rise to a three-peak structure in the quark spectrum 
just above T^. In Ref. 7), the Nambu-Jona-Lasinio model is employed as a chiral model, 
and the quark spectral function is calculated incorporating the effect of the dynamically 
generated soft modes into the quark propagator. It is noteworthy that although the soft 
mode obtained in the dynamical model employed in Ref. 7) is a composite particle with a 
finite width and thus not an elementary boson, a clear three-peak structure is realized in the 
quark spectral function. 

It would be interesting to investigate the effect in the change of the quark spectrum on 
various observables in the QGP phase. In particular, it would be intriguing to explore the 
dilepton production rate; because the quasiparticle peaks in the three-peak structure have 
low sound velocities, there can appear a sharp peak in the dilepton production rate owing 
to a mechanism similar to that causing the van Hove singularities. ^'^^ It would also be an 
interesting to consider the properties of mesonic excitations in the QGP phase composed of 
quarks having the spectrum obtained in this work. 

We remark that a lattice QCD simulation of the quark spectrum would be, of course. 



30 



preferable, although this would require the performance of several tasks, such as the gauge 
fixing, carrying out the chiral extrapolation and using dynamical (chiral) fermions to eluci- 
date the effects of the chiral soft modes. 

Contributions of higher order terms such as more than one-loop diagrams, which are 
not taken into account in the present, may affect the peak structure of the spectrum. In 
condensed matter physics, some corrections based on self-consistency conditions are pro- 
posed^*^) and their possible effects on the quark spectrum near the color superconducting 
phase transition are commented in Ref. 17). 

Because the present investigation is based on simple but generic Yukawa models with a 
massive scalar and vector boson, the results obtained here should apply to various systems 
at finite T composed of a light fermion and a massive boson, in addition to QCD matter 
in the QGP phase. One such system is the excitation spectra of neutrinos coupled to weak 
bosons at high T near the masses of the weak bosons. The investigation of a such system 
would be carried out exclusively in terms of the neutrino quasi-dispersion relations. 

In this paper, we have considered only massless quarks. A finite quark mass should affect 
the formation of the three-peak structure in the quark spectral function. Thus, the incor- 
poration of the quark mass effect is important for describing the chiral transition precisely, 
because the order of the transition, and hence, the quark spectrum near the critical point 
are sensitive to the quark mass. A finite quark mass also leads to a difference between the 
scalar (vector) and pseudoscalar (axial-vector) boson cases, and may suppress the three-peak 
structure, because it does not appear in Ref. 21), where a massive quark spectrum coupled 
with a massless boson is studied at finite T. Detailed study of such a effect is now under 
way and will be reported elsewhere. The finite chemical potential may also alter the 
three-peak structure because antiquarks in medium are fewer in this case, and thus, mixing 
between the quark and the antiquark-hole will be suppressed. Such a study is also now under 
investigation. 
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Appendix A 

Quark Self-Energy in the Yukawa Model with Scalar Coupling 

In this appendix, we summarize the calculation of the quark self-energy in the model 
with the scalar boson given in Eq. (3-2), 

/(Pk 

-—^gQ{k,iUn)V{p-k,iuJrn-i^n), (A-1) 
■II, ^ ' 

where 

Goik.iujn) = -. ^^-^ = J^Mkh'-. ^^, (A-2) 

[ii^n) - k -m^ ^ 2Eb{ii'n - tEb) 

are the free propagators of the quark and the scalar boson, respectively. Here, A±{k) — 
(1 ±7°7 ■ k)/2 are the projection operators onto positive and negative energies, respectively. 
Substituting Eqs. (A-2) and (A-3) into Eq. (A-1), we have 

d^k tA,{k)-f^ 1 1 



E(p,iu;m)^-g^J2^Yl J 

s,t=± n 



(27r)2 2Eb iujn — sEf iujm — i^n — tE^'' 

(A-4) 

with Ef — \k\ and E^ — y^(p — fc)^ -|- m^. Carrying out the summation over the Matsubara 
modes n and applying the analytic continuation icum cu -\-ir] in Eq. (A-4), we obtain the 
retarded self-energy and its imaginary part: 

u.)- V /■ ^vl.W/ f{sEf)+n{-tEb) 

E (p,a;)--5 2. / (27r)3 2E, u - sE. - tE, + irj^ ^^^^ 



s,t=± 



ImE^{p,u) = ng' J -0JA^ {f{sE,) + n{-tE,)] 6{u - sEf - tE,). 

(A-6) 



Equations (3-3) and (3-6) are obtained by carrying out the summations over s and t in 
Eqs. (A-5) and (A-6), respectively. 

We next decompose the self-energy with the projection operators into ^^(p, oj) = E+{p, uj)A^{p)'y^-\- 
Z'_(p, a;)yl_(p)7°. The imaginary part of U+{p,u!) is calculated in the following way: 

Imi;+(p,a;) = ilV[Imi:«(p, u;)^+(p)7°] 
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ng 



E 



s,t=± 



d?k 1 — sp ■ k 

(27r)3 AtEb 



{f{sEf) + n{-tEb)} 5{u; - sEf - tE^) 



(A-7) 



9^ 



-^Y^st dEf dE, [{\p\ - sEff - El + m"] 



^ s,t=± 

X {/(s^/) + n{-tE^)] 5(a; - sE^ - tE^) 



/■e+ 

X -a;)} / dEb5(u; - sEf - tE^), (A-8) 

^ e_ 



where in the third equahty we have used the relations 



p ■ k 

d^k 



p^ + E'j-El + m'' 



(27r)3 47r2 



2\p\Ef 
l^^j dEfEjJ^ ^ dE,E, 



(A-9) 
(A-10) 



with e± = \J {\p\ ± EfY + w?. 

In the time-hke region, one of the four terms in Eq. (A-8) takes a finite value (see Fig. 2). 
After some algebra, it is found that all these terms in the time-like region are reduced to the 
following simple form: 

Imr+(p, u) = Jj dEf [{\p\ +u- 2Ef){\p\ m^] 



with 



x{f{sEf)+n{sEf-u)}, 



111 

Ef = 



(A-11) 



(A-12) 



2(a;±|p|) 

In the space-like region, the two terms corresponding to (II) and (III) in Fig. 3 take finite 
values, and Eq. (A-8) yields 



ImZ".!. (p, uj) 



327rp2 



p — oo p1 

dEf- 

JEt Joo 



dEf 



x{f{Ef)+n{Ef-u)} 



327rp2 



El 



dEf - / dE 



[{\p\+uj-2Ef){\p\-uj) + m'] 



[{\p\+u;-2Ef){\p\-u;) + m^] 
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{f{Ef)+n{Ef-u;)} 





r^7 


327rp2 






\P\ - 
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p2 



327r p2 
where we have used the integral formulae 



(A-13) 



J — c 



de[f{e)+n{€-u;)]^ -co, (A-14) 

j ^ dee [f{e) - n{e - c.)] = ^ + -. (A-15) 

Equations (A- 11) and (A-13) are combined to give 
ImZ'+(p, cu) 

= £^ dE, [(Ipl +w- 2B;)(|p| - a,) + m=] {/(£,) + - u-)} 

For p = 0, from Eq. (A-7) we obtain 

Imi7± 0,u;) = coth— — + tanh — — . A-17) 

The T- independent part of ImZ'+(p, cu) is given by the T — hmit of Eq. (A- 16). In this 
case, the integral in Eq. (A- 16) can be performed analytically, and we obtain 

fl2 (p2_ 2\2 

lmE+ip,u;)T=o = - |p|) ^ ^4 sgnMg(P^ - m% (A-18) 

with P2 ^ - 

The T-dependent part, lmE^{p,ijj)T^o, is obtained through the replacement of the dis- 
tribution function in the curly brackets in Eq. (A-8) as 

f{Ef) + n{u; - Ef) n{u; - Ef) + f{Ef) + e{u;)e{P^ - rri^). (A-19) 

Appendix B 

Quark Self-Energy in the Yukawa Model with Vector Coupling 



The self-energy of the quark in a Yukawa model with a massive vector boson is given by 
Eq. (4-3) as 

E{p,iUm) = -/^X/ j J^:^^^^o{k,iuJn)Y'^i,Ap - k,iUm - iuJn): (B-1) 
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with the Proca propagator defined in Eq. (4-4). Carrying out the summation over the 
Matsubara modes n and applying the analytic continuation iuj„i uj + if], the retarded 
self-energy E^{p,uj) is obtained as Eq. (4-5), and the imaginary part is given by 

ImE (p,u) = -ng J^J j J^s^E, " ^) 

X {fisEf) + n{-tE,)} d{uj - sEf - tE,), (B-2) 

with Ef — \k\, El) — •y/(p — + and = {tEb,p — k). The quark sector of Eq. (B-2) 
gives 

hnE+(p,u,) = Tr [lmr«(p, a-)^+(p)7°] 



„ t L J 



dPk t 



{tEb — sk ■ q) {tEi, — p q) 1 — sk ■ p 



s,t=± 

X {fisEf) + ni-tEb)} 5iu - sEf - tE^) 



32Trp'^m 



' ^ s,t=± 



/■e+ 

+2sEf{p'^ - a;^ + 2m^)(|p| - uj)] / dEbS{uj - sEf - tE^) 

J e~ 



9 



2 rE 



327rp2m2 J ^ 



de[-{P'' - m'){{\p\ - ujf - 2m'} 



+2e{uj - \p\){P' - 2m')] {f{e) + n{e - cu)} 
;9{-P')l{P' - 2m'){u; - \p])n'T' 



327rp2777,2 

+u;[2m^ - P'{\p\{u; - \p\) + m'}] }. (B-3) 



where in the third equality we have used Eqs. (A-9), (A- 10) and 

. p'-Ej-Ei + m' ^ p'-E} + El-m' 

-2E, ' ^-^^ -2^\ ■ (^-'^ 

For p = 0, we have 

^ ^ ,^ s q' uj' + 2m' (u' — m')' ( ^ u' + m' , lj' — m' 
ImE±{0,u) = -^TTz^ 2 ^ 4 — ( ^ 



647r m2 cu^ \ 4Tuj 4Tuj 

(B-5) 

For the T-indcpendent part, which is obtained by taking the T — >■ limit in Eq. (B-3), 

we can perform the integral in Eq. (B-3) analytically, and we obtain 

a' P' + 2m' ( P' - m')' 

Imr^(p,c.)^.o = -|^(a; - IpI) ^ ^ gn(a;)g(P^ - m'), (B-6) 

with P'^uj'-p'. 
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Appendix C 

Renormalization of the T -Independent Part 



In this appendix, we perform the renormahzation of the T-independent part of the self- 
energy, E{p,u!)t=o, in the models studied in §§3 and 4. Because the imaginary part of 
E^{p,u!)t=o in both models is free from ultraviolet divergence, as shown in Appendices 
A and B, renormalization is needed for the real part only. The T-independent part is 
Lorentz covariant. To calculate S{p, (^)t=o by using this symmetry explicitly, we employ 
the self-energy in the Feynman propagator instead of the retarded one, L^. Once the 
renormalized self-energy in the Feynman propagator is derived, the corresponding retarded 
one is readily obtained. In the following, we omit the subscript of the self-energy, writing 
^Zy^-i—Q as ^ 

Because of the Lorentz invariance, the self-energy has the form 

E^(p,u;) ^i>E({P') + (P2), (C-1) 

with = {ui,p). For the self-energy used in §3 and §4, Z'2(P^) vanishes. 

We regularize the real part of the self-energy by employing the subtracted dispersion 
relation. This is much more convenient than other regularization procedures, e.g., the Pauli- 
Villars and dimensional regularizations. (We have confirmed that these regularizations give 
the same result in our models.) The subtracted dispersion relation is written 

where a denotes the subtraction point or the renormalization point, and Sthr denotes the 
threshold, which is given by Sthr = in our case. The parameters q are the subtraction 
constants, which are to be determined by the renormalization condition. Equation (C-2) is 
valid for a < s. 

C.l. Renormalization in the Yukawa model with scalar coupling 
From Appendix A, the imaginary part of U[ (P^) reads 



III 



2\2 



Imrf (P^) = p, ' 0{P' - m% (C-3) 

For the real part of Z'f , we use the subtracted dispersion relation given in Eq. (C-2) for n = 1 
to regularize the ultraviolet divergence. We impose the on-shell renormalization condition, 
i:f (p2 = a) = and a = 0. From Eq. (C-3), we have Im^f (P^ = 0) = 0. The renormalized 
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real part is then given by 



p2 

Rerf (P^) = — P 

TT 



ds 



Imrf(s) 
s{s-P^) 



327r2 



2\2 



log 



p2_^2 



3 



(C-4) 



The retarded self-energy, I]^[p,u!), is given by IieS^{p,u!) — IieIJ^{P'^) and 
Imi:^(p,u;) = e(u;)Imi:^(p2): 

.2 r (p2 _ ,„2^2 



r«(p,a;) 



log 



9 

32n^^ \ P4 

-^—P m e{u)e{P -m). 



p2 _ „^2 


3 






~ 2 ' 





(C-5) 



327r'' P4 

C.2. Renormalization in the Yukawa model with vector coupling 
Prom Eq. (B-6), the imaginary part of the self-energy E[ reads 

To renormalize the real part, we need the subtracted dispersion relation appearing in Eq. (C-2) 
for n = 2, 



p4 rco 

Rei:f(p2) = Co + CiP2 + — P / ds 



Imrf(g) 

S2(s-P2) 



Co + CiP^ 



5 



327r2 



2\2 



P2 + 2m' (P2 - m^) 



p4 



log 



p2_^2 



5 P2 ^ 
2 + 2 — 

6 777.2 p2 



(C-7) 



where we have set a = 0. We finds that the on-shcll renormalization condition Z'f(p2 
0) = is not sufficient to remove all the divergent parts. Here we impose an additional 
condition, ci = 0, to obtain the renormalized self-energy. 



Kerf (P2) 



9' 



327r2 



2\2 



P2 + 2m^ (P2 - m^) 



777^ 



p4 



log 



P^-m^ 



777^ 



5P2 ^ ^77^2 

2 2 — 

6 777^ P2 



(C-8) 



Therefore, the renormalized retarded self-energy is given by 

■log 



rf (p2) - ^ J ^ 



327r2 



771^ 



p4 



p2_^2 



771^ 



5P2 7^2 

2 + 2 — 

6 7772 p2 



. ^2 p2 + 277^2 (P2 - ^2)2 

^ ^ :r:; ^Sgn(a;)^(P'' - 777^'). 



327r 7772 



p4 



(C-9) 
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Appendix D 

Limiting Behavior of U^ijp.u) in the Yukawa Model 
with Vector Coupling 



In §4, we showed that the thermal mass of the quark coupled with the vector boson in 
the Proca formalism approaches uj — -\/6m in the high T limit, where m is the mass of the 
vector boson. In this section, we prove this limiting behavior. Because we are interested in 
the thermal mass, we only consider the quark propagator at vanishing momentum and omit 
the subscript and the argument of the momentum in the self-energy, writing E±{p — 0, a;) 
as in the following. 

Let us first derive the analytic form of Y{.eE{uj) in the high T limit, where the T- 
independent part is negligible. Thus, we consider the T-dependent part KeE{ijj)T^Q, which 
is calculated from the unsubtracted dispersion relation in Eq. (3-5), 



1 

Rer(a;)T^o = P / 



dz- 



cu — z 



uj^ — z^ 



(D-1) 



In the second equality, we have used the relation \mE(uj) = lmE{—uj), which is valid for a 
massless quark of vanishing momentum. Here, lmE{uj)T^o is given by 

g"^ ij? + Iw? [u? — m'^y 



647r 



X ( coth 

V 4Ta; 



+ tanh 



— 



m 



Substituting Eq. (D-2) into Eq. (D-1), we have the following: 



(D-2) 



327r2 



(/1H + 72H), 



(D-3) 



I'm 

h{u) = P / dz 
Jo 



1 z^ + 2m^ {z^ - m^^^ 
cu'^ — z^ m? 



coth 



z"^ + 
ATz 



tV /-i/" dx (2 + aV)(l-aV)2 



.(^) ^ P / 

J rr 



UJ ) Jq 1 — {af3xY 

X 

\Ax 

1 + 2m^ (z^ - ra^Y 



, / 1 a'^x\ , / 1 a'^x 

coth \ — - tanh — 

V4a; 4 / V 4a; 4 



dz 



up' — z'^ rn^ 



coth 



T 



2 r-Xja 



z^ + vr? 
ATz 

.,2„,2\2 



dy (l + 2aV)(l-«V) 



m J Jq 1 — a^'ip' I fP 

X 



coth + — ^ + tanh ;^ 1 - 2 

Ay A J \Ay A 



tanh 



z^ - m"^ 
ATz 



(D.4) 



, 2;^ — 

tanh 2 

ATz 



(D.5) 
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with a = m/T, (3 = m/uj, x = z/{ma) and y = m/ (az). 

In the high T hmit (a — > 0), the asymptotic forms of Ii and I2 are given by 




We have confirmed that the correction terms due to the finite value of a are 0(Q;ln(l/Q;)). 
Therefore, the self-energy in the high T limit is given by 

lim ReE(u;) = ^ (-- ^] . (D-8) 

In this case, the quasi-dispersion relation is given by the solution of the equation cu — 
ReE{u;) ^ —ReE{u;) — 0, i.e., the zero of ReE{u;). Therefore, in this model, the energy 
uj — y/6m gives the thermal mass in the high T hmit. 
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